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ABSTRACT 

Context. The highest redshift quasars at z ^ 6 receive considerable attention since they provide strong constraints on the growth of 
the earliest supermassive black holes. They also probe the epoch of reionisation and serve as "lighthouses" to illuminate the space 
between them and the observer. The source J1427+3312 (z = 6.12) has recently been identified as the first and so far the only known 
radio-loud quasar at z > 6. 

Aims. We investigated the compact radio structure of J1427+3312 on milli-arcsecond (mas) angular scales, to compare it with that of 
the second most distant radio-loud quasar J0836+0054 (z = 5.77) and with lower-redshift radio quasars in general. 
Methods. We observed J1427+3312 in phase-reference mode with ten antennas of the European Very Long Baseline Interferometry 
(VLBI) Network (EVN) at 1 .6 GHz on 1 1 March 2007 and at 5 GHz on 3 March 2007. 

Results. The source was clearly detected at both frequencies. At 1 .6 GHz, it shows a prominent double structure. The two components 
are separated by 28.3 mas, corresponding to a projected linear distance of ~ 160 pc. Both components with sub-mjy flux densities 
appear resolved. In the position of the brightest component at 1.6 GHz, we detected mas-scale radio emission at 5 GHz as well. The 
radio spectrum of this feature is steep. The double structure and the separation of the components of J1427+3312 are similar to those 
of the young (S 10 4 yr) compact symmetric objects (CSOs). 
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1. Introduction 

That the highest redshift quasars actually exist is their main in- 
terest. They provide significant constraints on the growth of the 
earliest supermassive (~ 10 9 M Q ) black holes. They probe the 
Universe close to the epoch of reionisation, and serve as "light- 
houses" that illuminate the space between them and the observer. 
Until recently (e.g. Springel et al. |2005l l, structure formation 
models in the ACDM Universe have been unable to produce 
massive quasars and galaxies whose the nuclei became active 
well before 1 Gyr following the Big Bang. Due to the small num- 
ber [just over 20, Jiang et al. (2008) and references therein] of 
z ~ 6 quasars known (of which only two are radio-loud), it is 
still an open question whether their physical properties are in 
general similar to those of lower redshift quasars. At present, we 
know little about quasars at z ~ 6, especially about the radio- 
loud ones, simply because our sample is very small. There are 
arguments (e.g. Haiman et al. 2004) in favour of a high surface 
density for z > 6 radio-loud quasars. 

The quasar FIRST J1427385+331 241 (J 1427+3312 here- 
after) was identified by McGreer et al. (2006) as the first radio- 
loud quasar at a redshift of z > 6. An independent identifica- 
tion was made by Stern et al. (2007 ). The source is found in the 
NOAO Deep Wide-Field Survey (NDWFS) region. Earlier ra- 
dio measurements by the Very Large Array (VLA) in the Faint 
Images of the Radio Sky at Twenty-centimeters (FIRST) sur- 



vey (White et al. I1997K the European Large- Area ISO Survey 
(ELAIS; Ciliegi et al. 1 19991 ). and by t he We sterbork Synthesis 
Radio Telescope (WSRT; de Vries et al. 120021 ) gave 1 .4-GHz flux 
densities of 1.7 - 2.1 mJy for J1427+3312. The source has been 
detected in several infrared bands. The optical and near-infrared 
colours, the broad absorption lines, and the possible intrinsic X- 
ray absorption implied by the Chandra non-detection (Murray et 
al. 120051 ) all suggest that dust extinction plays an important role 
in determining the appearance of this distant quasar (McGreer 
et al. 120061 ). Becker et al. ( 120001 ) find that the majority of their 
sample of radio-emitting broad absorption line (BAL) quasars 
selected from the FIRST survey are compact on the (X'2 scale. 
Apparently the BAL quasars are not like the normal quasars 
seen edge-on. The existence of a comparable number of flat- and 
steep-spectrum sources in the sample cannot be reconciled with 
simple unified models. Their preferred interpretation is that these 
objects are young or recently refueled quasars. 

The high-resolution radio structure of the second most dis- 
tant known radio-loud quasar (J0836+0054 at z — 5.77) is char- 
acterised by a single compact but somewhat resolved component 
having a steep radio spectrum between 1.6 and 5 GHz. Dual- 
frequency Very Long Baseline Interferometry (VLBI) observa- 
tions (Frey et al. 120031 120051 ) indicate that its radio emission 
originates from a region within ~ 40 pc. The aim of the VLBI 
observations presented here was to reveal the radio structure of 
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J1427+3312 on the linear scales of ~ 10 - 100 pc. Assuming a 
flat cosmological model with Ho = 70 km s _1 Mpc~', £2 m = 0.3, 
and = 0.7, the redshift z = 6.12 corresponds to 0.9 Gyr after 
the Big Bang (< 7% of the present age of the Universe). In this 
model, 1 mas angular separation is equivalent to a linear separa- 
tion of 5.65 pc. 



2. VLBI observations and data reduction 

We observed J1427+3312 with the European VLBI Network 
(EVN) at 1.6 GHz on 1 1 March 2007 and at 5 GHz on 3 March 
2007. The observed frequencies correspond to ~ 1 1 GHz and 
~ 36 GHz in the rest frame of the quasar. The a-priori coor- 
dinates with nominal uncertainties of 520 mas and 390 mas in 
right ascension and declination, respectively, were taken from 
Ciliegi et al. dl9991 l. At a recording rate of 1024 Mbit s -1 , ten 
antennas of the EVN participated in the 7-hour observations: 
Effelsberg (Germany), Hartebeesthoek (South Africa), Jodrell 
Bank Mk2 (UK), Medicina, Noto (Italy), Torun (Poland), Onsala 
(Sweden), Sheshan, Nanshan (P.R. China), and the phased array 
of the 14-element WSRT (The Netherlands). The data from the 
latter array were analysed as well, providing us with simultane- 
ous measurements of the total flux density of the source at both 
frequencies. Eight intermediate frequency channels (IFs) were 
used in both left and right circular polarisation. The total band- 
width was 128 MHz in both polarisations. The correlation of the 
recorded VLBI data took place later at the EVN Data Processor 
at the Joint Institute for VLBI in Europe (JIVE), Dwingeloo, the 
Netherlands. 

A mJy-level weak radio source, J1427+3312 was observed 
in phase-reference mode. This allowed us to increase the co- 
herent integration time spent on the target source and thus to 
improve the sensitivity of the observations. Phase-referencing 
involves regularly interleaving observations between the target 
source and a nearby, bright, and compact reference source (e.g. 
Beasley & Conway 1995). The delay, delay rate, and phase so- 
lutions derived for the phase-reference calibrator were interpo- 
lated and applied to J1427+3312 within the target-reference 
cycle time of 7 minutes. The target source was observed for 
~ 4.5-minute intervals in each cycle. The total observing time 
on J1427+3312 was 3.2 h at both frequencies. 

The phase-reference calibrator source we used was 
J 1422+3223, a compact extragalactic radio source with 1°36 
angular separation from J1427+3312. The J2000 right as- 
cension (a = 14 h 22 m 30 s .378956) and declination (6 = 
+32°23'10'.'44012) of the reference source in the International 
Celestial Reference Frame (ICRF) are available from the Very 
Long Baseline Array (VLBA) Calibrator SurvejQ. The posi- 
tional uncertainty is 0.49 mas. 

The US National Radio Astronomy Observatory (NRAO) 
Astronomical Image Processing System (AIPS, e.g. Diamond 
1 19951 1 was used for the data calibration and imaging. The vis- 
ibility amplitudes were calibrated using system temperatures 
regularly measured at the antennas. Fringe-fitting was per- 
formed for the calibrator and fringe-finder sources (J 1422+3223, 
J1407+2827, and J1331+3030) using 3-min solution inter- 
vals. The data were exported to the Caltech Difmap package 
(Shepherd et al. |1994| l for imaging. The conventional hybrid 
mapping procedure, involving several iterations of CLEANing 
and phase (then amplitude) self-calibration, resulted in the im- 
ages and brightness distribution models for the calibrators. 



1 http : //www . vlba . nrao . edu/ astro/calib/ index . shtml 
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Fig. 1. The naturally weighted 1.6-GHz VLBI image of 
J1427+3312. The positive contour levels increase by a factor of 
V2. The first contours are drawn at -50 and 50 /iJy/beam. The 
peak brightness is 460 /iJy/beam. The Gaussian restoring beam 
is 6.2 mas x 5.0 mas at a major axis position angle PA=29°. 
The coordinates are related to the brightness peak of which the 
absolute position is given in the text. 

Overall antenna gain correction factors (~ 15% or less) were de- 
termined and applied to the visibility amplitudes in AIPS. Then 
fringe-fitting was repeated in AIPS, now taking the clean com- 
ponent models of J1422+3223 into account. The residual phase 
corrections that resulted from the non-pointlike structure of the 
phase-reference calibrator were considered this way. The solu- 
tions obtained were interpolated and applied to the target source 
data. The visibility data of J1427+3312, unaveraged in time and 
frequency, were also exported to Difmap for imaging. The nat- 
urally weighted images at 1.6 GHz (Fig.Q]) and 5 GHz (Fig. [2]) 
were made after several cycles of CLEANing in Difmap. The 
lowest contours are drawn at ~ 3<x image noise levels. The the- 
oretical thermal noise values were 10 and 12 /^Jy/beam (lcr) at 
1 .6 GHz and 5 GHz, respectively. 

3. Results and discussion 

The source was clearly detected with the EVN at 1 .6 GHz, show- 
ing a prominent double structure (Fig.Q]). The two components 
are separated by 28.3 mas, corresponding to a projected linear 
distance of ~ 160 pc. Circular Gaussian brightness distribution 
model fitting in Difmap indicates that both components are re- 
solved (Table Q~|). The brightest one has 0.92 mJy flux density 
and 5.8 mas angular size (full width at half maximum, FWHM), 
and the weaker and more extended one has 0.62 mJy flux den- 
sity and 18.6 mas angular size. In the position of the brightest 
component at 1 .6 GHz, we detected mas-scale radio emission at 
5 GHz as well (Fig. [2j, fitted by a circular Gaussian model with 
0.46 mJy flux density and 3.2 mas angular size. The equatorial 
coordinates of the brightness peak are aj2ooo = 14 h 27 m 38'!58563 
and 5j2ooo = +33°12'41'.'9252. The positional uncertainty is de- 
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Fig. 2. The naturally weighted 5 -GHz VLBI image of 
J1427+3312. The positive contour levels increase by a factor 
of V2. The first contours are drawn at -50 and 50 /iJy/beam. 
The peak brightness is 167 yuJy/beam. The restoring beam is 
2.5 mas x 1.6 mas atPA=10°. 

Table 1. Fitted circular Gaussian brightness distribution model 
parameters of J1427+3312at 1.6 GHz (top) and 5 GHz (bottom), 
and estimates of the rest-frame brightness temperatures. 



Flux density 


Position (mas) 


Size FWHM 


Tb 


(mJy) 


North East 


(mas) 


(10 6 K) 


0.92 


1.37 -0.13 


5.75 


94 


0.62 


-24.73 -11.06 


18.62 


6 


0.46 


-0.29 0.22 


3.22 


15 



termined by the phase-reference calibrator source position accu- 
racy, the target-calibrator angular separation, the angular resolu- 
tion of the interferometer array, and the signal-to-noise ratio. In 
our case it is estimated as better than 1 mas, due to the accurate 
relative astrometry with respect to the reference source position 
provided by the technique of phase-referencing. 

The WSRT synthesis array data that were obtained during 
the VLBI observations were retrieved and processed with soft- 
ware tools developed at JIVE for the EXPReS e-VLBI project. 
Comparison of the total CLEAN flux density and the WSRT flux 
density measurements indicates that practically the entire radio 
emission of J1427+3312 originates from the components seen 
in our VLBI images. Based on the fitted flux densities, the radio 
spectral index of the brightest feature detected at both frequen- 
cies is a = -0.6 (S oc v", where S is the flux density and v is the 
frequency). J1427+3312 has also been observed and detected at 
1.4 G Hz with the VLB A in June 2007 (Momjian & Carilh |2007l 
CarilliUDnS}. The 1.6-GHz EVN image (Fig. [D is fully consis- 
tent with what is seen at 1.4 GHz with the VLBA (Carilli|2008). 
The significance of our dual-frequency EVN imaging results is 
that we can exclude the presence of a mas-scale flat-spectrum 
component in the source at ~ 50 /J.]y level. 

The radio-loud quasar with the second highest redshift 
known (J0836+0054, z = 5.77) has also been studied earlier 
with the EVN at 1.6 and 5 GHz (Frey et al. 120031 |2005l It 
shows a single compact component. Quite remarkably, its radio 
spectrum is also steep (a = -0.8). One could speculate that the 
two radio-loud z ~ 6 quasars known to date appear to resemble 



the well-known gigahertz peaked-spectrum (GPS) and the some- 
what larger compact steep-spectrum (CSS) sources (e.g. O'Dea 
1998). Indeed, the known parts of their radio spectra are simi- 
larly steep. A convex radio spectrum of J 1427+33 12 that could 
possibly be obtained with multi-frequency flux density measure- 
ments at lower frequencies would prove this assumption. The 
GPS/CSS sources are believed to be young (i.e. in the early stage 
of their evolution), or/and perhaps "frustrated" (i.e. confined by 
the dense interstellar medium). In the case of high-redshift radio 
galaxies, Blundell & Rawlings ( 1999) argue that the sources we 
see are necessarily young (< 10 7 yr) in a flux density-limited 
sample. This youth-redshift degeneracy is a consequence of the 
Malmquist bias (i.e. the higher the distance, the more luminous 
sources enter the sample) and of how the luminosity decreases 
with the source age. 

It is not surprising that both z ~ 6 radio quasars known to 
date have steep radio spectra. Early studies of the high-redshift 
(z > 3.5) quasar population predicted that the space density of 
flat-spectrum Doppler-boosted sources declines rapidly at high 
redshift (Savage & Peterson 1983). Investigating the evolution of 
the quasar radio luminosity function, Dunlop & Peacock (1990) 
showed an abrupt (a factor of ~ 5) cut-off of the comoving 
number density of the flat-spectrum radio quasars in the redshift 
range 2 < z < 4. A similar, albeit less pronounced, trend was 
found for steep-spectrum sources as well. Shaver et al. ((1996) 
argues for an even sharper decrease in the space density of radio- 
loud quasars at z > 3. More recently Jarvis & Rawlings (2000) 
showed that the sample of the most luminous flat-spectrum ob- 
jects is a mixture of Doppler-boosted and GPS sources. It com- 
plicates the interpretation of the redshift cut-off, which is proba- 
bly more gradual beyond z ~ 2.5. 

It is interesting to note that the double structure and the sep- 
aration of the components of J1427+3312 are similar to those of 
the young (< 10 4 yr) compact symmetric objects (CSOs), a sub- 
class of GPS sources typically found in radio galaxies at much 
lower redshifts of z S> 1 (e.g. Wilkinson et al. 1994, Polatidis 
et al. l2002l l. If J1427+3312 is indeed an unbeamed, young sym- 
metric object that we see in the early Universe, it might become 
possible to detect the expansion of its components with VLBI 
monitoring observations over the next decade or so. This could 
give us a straightforward tool for estimating the age of the ra- 
dio source. If we adopt the CSO model and a typical hot spot 
advance speed of ~ 0.3c derived by Taylor et al. (120001 ) for a 
sample of CSOs, we obtain a rough estimate of ~ 2000 yr for 
the kinematic age of the radio source in J1427+3312. 

The possibility that we see gravitationally lensed images of 
the same source at 1 .6 GHz cannot be completely ruled out by 
our observations. In fact McGreer et al. (2006) found two strong 
intervening Mg II absorption systems at z = 2.18 and z = 2.20 
in the spectrum of J1427+3312. Earlier simulation studies by 
e.g. Wyithe & Loeb (2002 ) suggest that a large fraction - up to 
one third, depending on the quasar luminosity function used - 
of the z ~ 6 quasars may be magnified by gravitational lensing 
by a factor of 10 or more. Despite the observational efforts, no 
case for strong gravitational lensing is found among the known 
z ~ 6 quasars so far. In the case of J1427+3312, a lens point 
mass of ~ 10 8 M located at z = 2.2 would be needed to pro- 
duce images that are separated by 28 mas. Strong Mg II ab- 
sorbers are known to be associated with "normal" field galax- 
ies (Steidel et al. |T994) that could well be harbouring a central 
black hole of this size. However, the impact parameters are gen- 
erally large, which suggests that the gravitational lensing effect 
of the Mg II systems may occur in the weak regime, enhanc- 
ing only the quasar magnitudes. Using a sample of nearly 7000 
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strong absorbers at 0.4 < z < 2.2, a recent study by Menard et 
al. ((2008 ) shows that even the gravitational magnification of the 
background quasars due to these systems is on average negligi- 
ble. We believe that strong gravitational lensing is unlikely to 
play a role in the case of J1427+3312. 
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4. Conclusions 

We imaged J1427+3312, the highest-redshift radio-loud quasar 
known to date with the EVN at two frequencies (1.6 and 5 GHz). 
At 1 .6 GHz, our phase-referenced VLBI image shows two re- 
solved radio components with peak brightness ratio ~ 3 : 1 and 
flux density ratio ~ 3 : 2. The brightest of the two components is 
also detected at 5 GHz, suggesting a steep radio spectrum with 
spectral index a = -0.6 between the two frequencies observed. 
We determined the accurate ICRF position of the source. The 
high-resolution radio structure of J1427+3312 is similar to that 
of the CSOs. This may imply a young age of the radio source, 
which could possibly be verified by future monitoring of the ex- 
pected change in the components' separation. High-sensitivity, 
high-resolution radio imaging of a larger sample of z ^ 6 radio- 
loud quasars still to be identified spectroscopically could provide 
information on how common these very compact steep-spectrum 
objects are in the early Universe. 
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